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We calculate the Casimir interaction between a sphere and a plate, both described by the plasma 
model, the Drude model, or generalizations of the two models. We compare the results at both zero 
and finite temperatures. At asymptotically large separations we obtain analytical results for the 
interaction that reveal a non-universal, i.e., material dependent interaction for the plasma model. 
The latter result contains the asymptotic interaction for Drude metals and perfect reflectors as 
different but universal limiting cases. This observation is related to the screening of a static magnetic 
field by a London superconductor. For small separations we find corrections to the proximity force 
approximation (PFA) that support correlations between geometry and material properties that are 
not captured by the Lifshitz theory. Our results at finite temperatures reveal for Drude metals a 
non-monotonic temperature dependence of the Casimir free energy and a negative entropy over a 
sizeable range of separations. 



I. INTRODUCTION 

The past decade has witnessed rapid progress in the 
precision of Casimir force measurements [1]-^] • The mea- 
surement precision that is expected in the near future 
demands accurate theoretical calculations of the Casimir 
force for the geometries and materials used in experi- 
ments. While Casimir's original calculation for perfect 
metal plates [7] and Lifshitz's formula for dielectric slabs 
[8] only apply to planar, parallel surfaces, recent measure- 
ments have set limits on geometry induced corrections in 
the most frequently used sphere-plate geometry [S]. The 
geometry dependence of Casimir forces is intriguing as 
it can vary substantially with the shape and relative ori- 
entation of the objects [TDHH]. Material dependence in 
the form of dissipation in metals has been experimentally 
confirmed to have an effect on the Casimir force [3 [13] . 
It is thus important that the geometry and material de- 
pendence be carefully investigated for the experimentally 
most important sphere-plate configuration. 

In order to compare the experimental results to theory, 
the Derjaguin or proximity force approximation (PFA) 
[H] has commonly been used. This approximation ne- 
glects the non-additivity of Casimir forces by estimating 
the interaction between curved surfaces in terms of the 
planar surface interaction between infinitesimal and par- 
allel surface elements. Its validity is hence limited to the 
singular limit of vanishingly small separations between 
the surfaces. A systematic extension to larger separa- 
tions is not possible within such approximations. 

The first exact computation of the Casimir interaction 
energy for a perfectly reflecting sphere and plate was pre- 
sented in Ref. [TS]. Recently, corrections that come from 
using the plasma or Drude model were computed at zero 
temperature [16] and at T = 300° [I?]. Other open ge- 
ometries with curvature such as a cylinder above a plate 



have been studied for perfect metals [T5] . Corrections to 
the PFA in the case of perfect metals for a cylinder above 
a plate and a sphere above a plate have been obtained us- 
ing path integral approaches [TH [20] and for scalar fields 
employing a world line formalism [5T]. 

Here we show that Casimir forces reveal a rich inter- 
play between geometry (radius of the sphere and ob- 
ject's separation), optical properties of metals and ther- 
mal fluctuations. We study this in detail by calculat- 
ing the Casimir interaction for different sphere radii and 
separations using the (i) the Drude model, (ii) a gen- 
eralized Drude model, (iii) the plasma model and (iv) a 
generalized plasma model at different temperatures. The 
study of these combined effects is of utmost importance 
since Casimir force measurements continue to be carried 
out using this geometry and an increasing accuracy is 
expected. Hence, the experimental findings will begin 
to show sensitivity to the material and temperature ef- 
fects, which we take into account here. Furthermore, the 
unabated controversy whether the plasma or the Drude 
model is more appropriate for describing the optical prop- 
erties of metals in Casimir calculations compels us to pro- 
vide results for both models so that experimentalists can 
build on them when studying this problem further. The 
plasma model is a high-frequency description of the opti- 
cal properties and the divergence ~ 1 /oj^ of its dielectric 
function for small w is unphysical for metals. The Drude 
model provides a proper low-frequency description for 
metals with a l/cj divergence of the dielectric function 
for small lo. At large frequencies, both models become 
identical. 

Below, we supply numerical results for the Casimir 
interaction at arbitrary separations as well as analytic 
formulas for the asymptotic interaction at large separa- 
tions. Depending on the model under consideration, the 
asymptotic results show universal or non-universal (i.e., 
material-independent or -dependent) behavior, a feature 
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which is not present for the simple case of two parallel 
metal plates and hence results from the interplay of finite 
object sizes and material properties. 

II. GENERAL EXPRESSION FOR THE 
INTERACTION 

To calculate the interaction of a metallic sphere of 
radius R and a metallic plate with a separation d be- 
tween the center of the sphere and the plate, we employ 
a scattering approach for Casimir interactions, which is 
described in detail in Ref. [TT]. The Casimir free energy 
of this system at temperature T is given by 

oo ' 

E = keTY, log det[l - F («:„)] , (1) 

with Matsubara wave numbers Kn = 2i:nkBT /he. The 
primed sum indicates that the contribution for n = is 



to be weighted by a factor of 1/2. At zero temperature 
the sum is replaced by an integral along the imaginary 
frequency axis [TT], 

he r°° 

E=— log det[l - Y{k)] , (2) 
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where the matrix Y is given by the product 

'-Iml'm' ~ s,lm^lml'm""i-m K"^) 

of the T-operator of the sphere and an operator U 
that describes the propagation of waves between the plate 
and the sphere and the scattering of them at the plate 
(see below). We represent these operators in a vector 
basis of spherical waves, where a, /3=E, M denote elec- 
tric or magnetic multipoles and I, m label the spherical 
waves. For a sphere of radius R with uniform permittiv- 
ity e{uj) and permeability ^(w) the T-matrix elements for 
M-multipoles are given by 
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with n — e{iK)iJ,{iK), rj — e{iK) / ii{iK) . The T- 
matrix elements for E-multipoles, T^i^^^, are obtained 
from Eq. Q by interchanging e and /i and by changing 
the overall sign. By taking e{iK) — > cx) at an arbitrarily 
fixed /i(iK) in Eq. Q, the limit of a perfectly reflecting 
sphere and plate is obtained. Then the matrix elements 
become independent of /i. 
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These matrix elements scale for small + . It is 

interesting to compare this behavior to the scaling of 
the general matrix elements of Eq. Q for the dielec- 
tric functions of the Drude and plasma model. For both 
models the Tf;^ ~ behavior is unchanged for E- 

multipoles. The coefficients become material (plasma fre- 
quency) dependent for the plasma model but retain the 
universal values of a prefect reflector for the Drude model. 
However, for M-multipoles only the plasma model shows 
this universal behavior while the Drude model yields a 
different scaling T^\^ ~ with non-universal con- 

ductivity dependent coefficients. 

The operator U can also be expressed in a spherical 
wave basis. It describes the propagation of waves from 



the sphere to the plate, a reflection at the plate and the 
propagation back to the sphere. The reflection of waves 
at a dielectric plane is described most easily in a plane 
wave basis with in-plane wave vector fcy . The T-matrix 
elements of the plane are then given by the usual Fres- 
nel coefficients. The conversion from plane to spherical 
waves and simultaneous translation from the sphere to 
the plane is obtained by multiplying the plane's T-matrix 
from left and right by a matrix £'hn(fc||). After defining 
z — fc|| /k, the matrix multiplication runs over the contin- 
uous variable z and the elements of the operator U can 
be written as 



U, 
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-J2Dr;:i{z)T^{n,z)D^,l*iz), 

7 

(6) 

where the plate's diagonal T-matrix, T^{k,z), for polar- 
ization 7 are given by 



_ _yu(iK)\/l + - y^e{iK)fi{iK) + z^ 

^ fl{iK)^/l + + y^e{iK)fJ,{iK) + z^ 



p ~ 



z(iK)\/l + — e{iK)n{iK) + z"^ 
z{iK)\/l + z^ + ^ t(iK)[i(iK) -I- 



(7a) 
(7b) 



The exponential factor in Eq. ([6]) describes the transla- 
tion from the sphere to the plane and back by a total 
distance 2c? in the plane wave basis. The elements of the 
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matrix that converts between plane and spherical waves 
are given by 



also for metals described by the plasma and Drude model 
at zero and finite temperatures. 
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III. LARGE DISTANCE INTERACTION AT 

T = 

In this section we consider the zero temperature 
Casimir interaction at large separations for different di- 
electric functions. 



where P™ is the associated Legendre polynomial of or- 
der /, m. These elements have the following symmetries 
under complex conjugation, 
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In what follows we employ Eqs. ([I]), ([2]) to obtain the 
Casimir interaction for perfectly reflecting bodies and 



Perfect reflector 



In the limit of perfect reflectivity of the plate, one finds 
from Eq. ([t]) with e — > cx) the simple T-matrix elements 
j*i _ rpE _ ^ With this simplification, the integration 
over z m Eq. ^ can be performed analytically. We find 
for the elements of U the same result that was obtained 
before, using the method of images [15] , 
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Using this result and the T-matrix elements of Eq. ^ 
we obtain for the interaction energy the large distance 
expansion 

9 i?3 
16 "d^ 



E = -- 



hc 



25 

32 d6" 



(12) 



at zero temperature |15j . 



B. Plasma model 

We now assume that both the sphere and the plate are 
described by the plasma model which on the imaginary 
frequency axis has the dielectric function 



epiicK) = 1 + 



27r 



(13) 



The plasma wavelength Ap is related to the plasma fre- 
quency ujp by Ap = 2'Kc/ujp. Note that the plasma model 



provides a high-frequency description of optical proper- 
ties and ignores dissipation. Hence it is not expected to 
capture the low frequency response of a metal. To under- 
stand the physical meaning of the results for the Casimir 
interaction presented below, it is interesting to realize 
that the dielectric function of Eq. ( 13 ) appears also in 



the wave equation for the magnetic field in a supercon- 
ductor when it is described by the London theory. The 
second London equation and the Maxwell equations yield 
tp with the penetration depth Ap — ^ rrieC^ / (IGn^nse"^) 
for superfluid carriers of density Us, charge e and mass 

TOg. 

To obtain the large distance behavior of the Casimir 
energy, we need to expand the T-matrices for small k. 
To this end, we set k — u/d and expand the relevant 
expressions in powers of 1/d. The T-matrix elements 
of the sphere scale as for k ^ for both E and 

M polarizations. In the case of the E polarization the 
coefficients are universal and are given by the perfect 
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reflector result which corresponds to 
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However, for the M polarization the coefficients are not 
universal and depend on the plasma wave length as fol- 
lows 



1 



/,_i(2^i?/Ap) (2/ + l)!!(2Z-l)!! 



(uR/d) 



21 + 1 



(15) 

In the limit of a small plasma wavelength, Xp <^ i?, the el- 
ements of this matrix approach the perfect reflector limit 
with is given by 
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For a large plasma wavelength, Ap ^ i?, the elements are 
not universal and reduced by a factor {R/Xpf' compared 
to the perfect reflector limit. 



(2/ + l)!!(2Z + 3)!! 



(17) 



The latter result can be understood in terms of the Lon- 
don superconductor interpretation of the plasma model. 



If the penetration depth Ap becomes much larger than 
the radius, the sphere becomes almost transparent for the 
magnetic field and the T-matrix elements are reduced to 
small values ~ i?^'+^/Ap. 

The T-matrix elements of the plate with ep{icK) of 
Eq. ( 13 ) depend also on the lateral wave vector kii . To 



obtain the large distance expansion, we set k\\ = v/d and 
expand the T-matrix for large d with z = fcy/K = v/u 
fixed. This yields the expansion of the plate's T-matrix 
elements. 
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With this expansion, the integral over z in Eq. ^ can 
be performed analytically, and one obtains an expansion 
in 1/d of the matrix elements of U which depend on u 
and Xp/d only. When we substitute the matrix elements 
of Eqs. ([U]), (|l5]) and ([l8]) into Eq. ^ and expand the 



energy of Eq. |TF in powers of 1 /d, we obtain the interac- 
tion to order ifd^ by including I 
result can be written as 



2 partial waves. The 



with the fimctions 
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(20a) 
(20b) 



(20c) 



Note the coefficient of the leading term depends on 
Ap and hence is not universal. Only in the two limits 
Ap/i? — )■ and XpjR — )■ oo the coefficient approaches 
the material independent values 9/16 and 3/8, respec- 
tively. This behavior is consistent with the two limiting 
forms of the sphere's T-matrix of Eqs. (16), ( [l7| . The 
limit XpjR — >■ describes perfect reflection of electric and 
magnetic fields at arbitrarily low frequencies and hence 



agrees with the result of Eq. ( 12 ) where for dipole flue 



tuations the E polarization yields twice the contribution 
of the M polarization, cf. Eq. (ImI and Eq. ([l6|) for Z 1. 



For XpjR — T' oo the coefficient /4 is reduced by a factor 
2/3 since the M polarization does not contribute to the 
leading term ~ R? jd'^ due its suppression by [R/XpY, 
cf. Eq. (17 1. Physically, the non-universal behavior of /4 



can be understood when the objects are considered as 
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London superconductors. For \p/R — >■ a static mag- 
netic field is perfectly screened and the objects become 
perfect reflectors. If Xp/R 1, a static magnetic field 
can penetrate the entire sphere and hence the M polar- 
ization does not contribute to the Casimir energy. From 
this interpretation it follows that normal metals, which 
can be penetrated by a static magnetic field, should in- 
teract to leading order in R/d only via E polarizations 
leading to fi = 3/8. We shall reach the same conclu- 
sion when we consider the Drude model below. The co- 
efficient of R'^/d^ is always positive and varies between 
(13/207r2)(Ap/i?) for \p/R and {3 /lOn^) (Xp/R) for 
Xp/R — > oo. The coefficient of R^/d^ can be negative 
(for Xp/R — > 0) or positive. In Sectioijvj we compare 
the exact findings of this Section to our results from a 
numerical evaluation of Eq. ([2| over a wide range of sep- 
arations. 

Finally, it is instructive to compare the above results 
to the interaction between two parallel and infinite plates 
that are described by the plasma model. In this case, 
the large distance expansion applies to d 3> Ap and the 
leading term is given by the universal perfect reflector re- 
sult. The plasma wavelength appears only in corrections 
to the leading term that can be expanded in powers of 
Xp/d. This universal behavior is a consequence of the 
(unrealistic) assumption of an infinite lateral size of the 
plates which removes any finite length scale of the object 
that could be compared to Xp. Hence, a finite penetra- 
tion depth only yields an increased effective separation 
which for d^ Xp obviously approaches d, explaining the 
universal large-d result. 



C. Drude model 

The Drude model describes the low-frequency response 
of a metal which depends on its dc conductivity a. 
For large frequencies it becomes identical to the plasma 
model with plasma wavelength Xp. On the imaginary 
frequency axis, the Drude dielectric function is given by 



eoiicn) = 1 + 



(2^)2 



[XpKY + TTCK/ap 



(21) 



The conductivity is associated with the length scale Xa = 
27rc/(T. At large distances d, we need to consider the limit 
of small K at fixed z — k]^\/K for the plate's T-matrix, 
which yields with k = u/ d 



P 



= 1 



Y irdap 2ndap ' \ J 

(22a) 

(22b) 
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T^dcTp ^ + 1 2TTd(7p 



The approach of unity for both polarizations is a conse- 
quence of keeping kw/n fixed in the limit k — ^ 0. This 



behavior arises from the fact that the plates are infinitely 
extended so that arbitrarily small fey are allowed. The sit- 
uation is different at finite temperatures where one has to 
take K — >■ at fixed fcy for the first term of the sum over 
Matsubara frequencies. In the latter limit the magnetic 
contribution T^ ^ vanishes. 

For the sphere with the Drude dielectric function of 
Eq. (21) we obtain for the T-matrix elements with I — 1 



the low frequency expansion 
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(23a) 
(23b) 



While the leading term of the E polarization agrees with 
the perfect refiector result, the leading term of the M 
polarization is reduced by a factor Rk — uR/d compared 
to the perfect reflector case. Therefore, one expects that 
only the E polarization contributes to the leading term 
of the interaction at large distances. 

With the above expansion of the T-matrix elements 
the integrations over u and z can be performed and from 
the dipole contributions with I = 1 we obtain for the 
energy the large distance expansion 
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The leading term in Eq. (24) shows the universal am- 



plitude coming only from the E polarization as expected 
from the form of the T-matrix elements. This result re- 
produces the prediction of the plasma model in the limit 
where Xp R, see the discussion below Eq. (20a). This 



limit describes the situation where a static magnetic field 
can fully penetrate the sphere and hence describes a nor- 
mal metal. The correlations between material and shape 
become obvious when one compares the above result to 
the interaction between two parallel and infinite plates 
that are described by the Drude model. For this geome- 
try the large distance expansion applies to d ^ c/a. The 
leading term of this expansion is identical to the prefect 
reflector result, as for the plasma model. The dc con- 
ductivity appears only in corrections to the leading term 
that can be expanded in integer powers of \fcfad. Since 
the frequently used PFA for the sphere-plate geometry is 
based on the two-plate energy, it would predict at suffl- 
ciently large d for both the plasma and the Drude model 
the perfect reflector result which has equal contributions 
from E and M polarization. However, it is known that 
the PFA does not apply to large distances. It should 
be noted that the result of Eq. ( 24 ) cannot be applied to 



an arbitrarily large dc conductivity a since then the term 
^ i?^, which comes from the M polarization of the sphere, 
diverges. The condition for the validity of Eq. (24) can 
be written as c? ^ i?, Ao-, R^ /X^,. Below we shall study 
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the vahdity range of this expansion further by comparing 
it to numerical results. 



IV. HIGH TEMPERATURE LIMIT 

In this section, we study the high temperature limit 
of the sphere-plate interaction for the plasma and Drude 
model. In this case, the interaction is given by the first 
term of the Matsubara sum of Eq. ([T]). Hence we have 
to compute the matrix elements of Y{ko = 0). This 
zero-frequency result will turn out to be also useful when 
computing the Casimir energy at zero and finite tempera- 
tures below since the limit k — >■ is numerically unstable 
due to the divergence of certain Bessel functions. 



A. Plasma model 

Here we have to consider the limit k — >■ at fixed fc|| of 
the T-matrix of the plate since we are interested in arbi- 
trary separations d. In this limit the T-matrix elements 
are given by 



rpM _ _ 



= 1 



(25) 



The elements for the M polarization are non-universal 
and vary between 1 for Ap — >■ (perfect reflector) and 
for Ap — >■ oo. The latter limit can be interpreted as a 
London superconductor with diverging penetration depth 
such that the plate is transparent to a static magnetic 
field. 



for 



K 



For Ap ^ the U-inatrix of Eq. ^ can be obtained 



from Eq. (10 1. To obtain the U-matrix for 
K — ?> at non-zero but small Ap <C d, we set fcy = v/d 



and expand T^f^ of Eq. (251 in Xp/d so that the integral 
of Eq. ([6| can be perforniea analytically. Since we are in- 
terested in the limit k — > 0, we only need the conversion 
matrix elements -D;"^(z) for large arguments z. At large 
z the associated Legendre poynomials Pim{z) assume the 
limiting form {—i)™'{2l — l)!!z'/(Z — to)!. Using the inte- 
gral e~^'"v"'dv = n!/2"+-'^ we obtain to leading order 
for small k the matrix elements 
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(26a) 
(26b) 



-i-i' 



where the matrix elements of Y are given by Eqs. (|3|, 
and hence can be ignored. The T-matrix ele- ( 14 1, ( 15 ) and ( 26 ). By truncating the matrix Y at lowest 



The matrix elements U^,^, U^^f,^ scale for small n as 



{ndy 

ments of the sphere for k — are given by Eqs. ( 14 ) , ( 15 ) 
and hence scale as k^'"*"^. The 1ow-k scaling of the ma- 
trix elements of U and Tg shows that the elements of 
the matrix Y scale as y^^f^ ~ Yi'^^l ~ k'"'' and 

Yil^L ~ yl'lv.^ - Hence for ^ the cou- 

pling of E and M polarization does not contribute to the 
energy. We set again k — u/d and introduce the rescaled 

matrix Y with elements Yimi'm = u^^YimVmU^ so that 
divergences for it — > are removed and in that limit all 
elements of and Y^^ assume non-zero finite val- 

ues that depend on R/d and Xp/d, and all elements of 
Y^^ and Y^^ vanish. The rescaling does not change 
the determinant of Eq. ([T]) so that det Y = det Y . In the 
high temperature limit the energy can then be written as 



hj = —— log det 



1 





f^^(w-^O) 
(27) 



order ^ = 1 we get the high-temperature free energy 
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27r 
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R^ 
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(28) 



which applies for d ^ R, Ap, Ay — hc/kBT. Notice that 
this energy is not universal in the sense that the leading 
term depends on the plasma wave length. For Ap 
R, the amplitude of the leading term becomes —3/8, in 
agreement with the high-temperature result for perfect 
reflectors [TT] . For Ap ^ R, the amplitude of the leading 
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term approaches —1/4 which is identical to the result 
for the Drude model (see Eq. (27 1 below). The behavior 



in these two limits is consistent with the corresponding 
limits of the zero-temperature result of Eq. (19). 



B. Drude model 

For this model, the T-matrix of the plane for k — > at 
fixed k\\ behaves differently from Eq. (22 1. While — 1, 
the magnetic part vanishes, T^' -> 0. Eq. ^ shows 
that to leading order for small k, the matrix elements 



^;mf'm are given by Eq. (pel) 



In fact, we do not need 
to find the other matrix elements of U : the elements 
coupling unlike polarizations are reduced by a factor k, 
and the elements U/!^ff^ 



are multiplied by Tj^^ of the 



sphere, which scales as k^'"*"^ for small values of k, and are 
thus smaller by a factor k also. The (universal) elements 
of Tfi^^ for small k are given by Eq. (14 1. This shows 



that only the E polarization contributes to the energy at 
high temperatures and from Eqs. (26) and (14 1 follows 



the explicit result for the elements of the rescaled matrix 



l + l 1 



ll{2l + l)l'{2l' + 1) {2V - + 



(R/d) 



21 + 1 
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(29) 



In the high temperature limit the energy is then given by Eq. ([T]) carefully for k — )■ 0. In Sect. IV we have already 



i?=M:iogdet[l 



(30) 



Notice that this result is universal at all separations since 
the matrix Y depends only on R/d. The absence of mag- 
netic contributions is in agreement with the high tem- 
perature interaction between two parallel plates that are 
described by the Drude model. A truncation of the ma- 
trix F at / = 2 and expansion of Eq. (30) for small R/d 
yields the large distance result 



E : 



-kuT 



1 /rV 1 fRV 3 



4:\d J 4:\d 
which applies when d^ R, Xt- 



R\ /I 
128 \ d) ^'^[dJ 



(31) 



V. NUMERICS 

In this section, we evaluate the Casimir energy based 
on Eq. ^ for zero temperature and Eq. ([T]) for finite tem- 
peratures. Our results are obtained by numerical com- 
putation of the determinant, the integral over k (or sum 
over n) and the integral over z of Eq. ([6]). The matrix 
Y is truncated at a finite partial wave order, £max- We 
chose Imax such that the result for the energy changes 
by less than a factor of 1.0001 upon increasing imax by 
10. The required value of imax depends on the separa- 
tion between the plate and sphere. As the separation 
decreases, imax increases. For example for R/d < 0.75, 
we used imax = 24, whereas for R/d — 0.8 and 0.85, one 
needs imax — 34 and for R/d = 0.9 the value imax = 54. 

The numerical computation of the determinant, the in- 
tegrals and sum poses no principle problem. However, it 
is important to consider the determinant of Eq. ^ or 



seen that the matrix elements Yu 



for small k scale as 



j-i' 



J-i'+i 



K' ' or K' " ' This shows that for small k the matrix 
elements with I 3> I' become extremely small whereas 
those with I ^ /' increase rapidly. For large values of 
Imax this behavior makes the computation of the deter- 
minant at K = numerically ill-conditioned. However, 
the analytical results presented in Sect. IV allow us to 
calculate the n = term in Eq. (jlj or the integrand of 
Eq. ([2]) at K = for imax — 100 and even larger. In fact, 
as the value of is increased beyond 0.9, larger values 
imax > 60 must be used in order to accurately calculate 
the energy. For sufficiently high temperatures, the sec- 
ond Matsubara wave vector ki = 2'KkBT /he in Eq. ([T]) 
becomes sufficiently large and hence poses no numerical 
problem for the computation of the energy. For exam- 
ple, for T = 300° iiT the Casimir energy can be calculated 
for R/d — 0.95 with imax — 72. As imax increases, the 
interval in the vicinity of k = in which the integrand 
cannot be obtained with sufficient precision numerically 
increases too. Due to this behavior, we restrict the cal- 
culation of the Casimir energy at T = to R/d < 0.9. 



A. Casimir interaction at T = 

In this section, we calculate the Casimir energies for 
the usual Drude and plasma model given in Eqs. ( 13 ) 



and (21), respectively, for parameters of gold as given 



below. There are three dimensionless parameters which 
we choose as d/R, Xp/R and Xp/Xa- The first two pa- 
rameters can be controlled for a given material by chang- 
ing the separation d and the radius R of the sphere. In 
order to avoid strong finite size effects in the electronic 
response, we assume that Ao 



Xp < R. 



In Sections IIIB and IIIC using i. 



par- 
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tial waves, we obtained an asymptotic expansion of the 
Casimir energy for both plasma and Drude model at 
large separations, see Eqs. (19) and (24). In Fig. [T] we 
compare the analytical results to the exact numerical re- 
sults that were obtained as described before. The graph 
shows the exact energies for the Drude and the plasma 
model normalized to the exact energies for perfect re- 
flecting surfaces, taken from Ref. [15]. For the plasma 
model we used \p/R = 0.05 and 0.5, respectively, and 
for the Drude model the same two values for \p/R and 
we set Xp/Xcr = 27. 4. The figure illustrates the material 
dependence of the Casimir energies. For large separa- 
tions, the ratios for the plasma model approach values 
slightly smaller than one, which is consistent with the 
Ap/i?-dependent asymptotic form predicted by Eq. (19 1. 
For the case of the Drude model, the ratio tends to the 
universal number 2/3 at large separations, as predicted 
by Eqs. (12), (24). In the case of the plasma model, the 
asymptotic result describes the energy up to R/d « 0.4 
nicely. For the Drude model, however, the agreement be- 
tween the analytical and numerical findings is limited to 
extremely small R/d < 10~^. This example clearly indi- 
cates distinct correlations between material and geome- 
try. Our result shows that for the Drude model a larger 
number of partial waves than for the plasma model is 
necessary to accurately calculate the Casimir energy. 




0.5 - 



0.4 - 



FIG. 1: i?/i?perf.refl. agaiust R/d for the plasma model for 
\p/R = 0.5, 0.05 (open circles and filled squares respetively) , 
and the Drude model for the same values oi \p/R (filled circles 
and open squares, respectively) and Xp/Xa ~ 27 A. The solid 
and dashed lines represent the asymptotic results of Eqs. ( 19 1, 



([24| for the plasma and the Drude model, respectively. Inset: 
Magnification of the short distance range. 



We also compare the exact numerical results with the 
Casimir energy obtained by the PEA for both the plasma 
and the Drude model. The PEA energy is obtained by 
integrating the PEA force F = 2TrREp\atcs{d— R) with re- 
spect to d, where Epig^tes{_d,—R) is the energy of two paral- 
lel plates at distance d—R as given by the Lifshitz formula 
[22] with the corresponding dielectric function. Fig. [2] 
shows the exact Casimir energy calculated numerically 
for the plasma model with plasma wavelength Xp = 0.05 
and Xp — 0.5, respectively. The figure shows also the 
PEA energy for the same values of Xp. As expected, the 
discrepancy between the exact and PEA energy decreases 
as R/d increases and is expected to vanish for d — )■ i?. 
This is clearly visible from Fig. [3] which shows the rela- 
tive corrections to the PEA energy at short separations. 
Interestingly, the dependence of the corrections on Xp is 
not fully described by the Lifshitz theory since the data 
for different Xp do not collapse onto a single curve. This 
demonstrates correlations between geometry and mate- 
rial properties that are not described by the PEA. For 
example, for Xp/R = 0.5 and Xp/R — 0.05 we find at the 
shortest studied separation oi d — R = O.lli? the exact 
energy to be 85% and 87% of the PEA energy, respec- 
tively. For perfect refiectors the reduction was found to 
be « 87% at the same distance [T3|. We find similar 
results using the Drude model. The energies associated 
with the Drude model are not shown here since they col- 
lapse onto the data for the plasma model at short sepa- 
rations. 




0.2 0.4 0.6 0.8 

R/d 



FIG. 2: Numerical result for ~Ed/hc against R/d for the 
plasma model at Xp/R — 0.05 (squares) and at Xp/R — 0.5 
(circles). The lines represent the PFA energy at Xp/R = 0.05 
(dashed) and at Xp/R = 0.5 (solid). 



B. Casimir interaction at T 7^ 

The Casimir free energy at finite temperatures T de- 
pends on the thermal wavelength At — hc/ksT. This 
additional length scale introduces an additional dimen- 
sionless parameter Xt/R. To investigate the infiucnce 
of temperature, we calculated the Casimir free energy 
at two different values of this parameter. We have cho- 
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0.4 0.6 O.i 

(d-R)/R 

FIG. 3: Ratio of the numerical results for the Casimir energy 
shown in Fig. [2] and the PFA energy based on the Lifshitz 
theory for the plasma model with Xp/R = 0.05 (squares) and 
Xp/R = 0.5 (circles). The ratio is shown as a function of the 
surface-to-surface distance d — R. 



are the resonant frequencies of the core electrons. Typ- 
ical parameters for gold are given by Wp = 9eV for the 
plasma frequency and 7 — 35meV for the relaxation rate 
|16) . These parameters correspond to the length scales 



Ap = 2-Kc/ujp = 137nm and 



— 27rc/7 = 35.4//m. For 
these parameters, the dc conductivity a — uj'^/{At:^) is 
184. 2eV, corresponding to the length scale A^. = 27rc/(T = 
6.7nm. Note that electron scattering is not described 
by the usual plasma model. However, as can be seen 
from Eq. ( 32 ) , dissipation is included in the generalized 



plasma model due to the interband transition of core 
electrons. To calculate the Casimir energy, we use the 
oscillator parameters of gold which are presented in Ta- 
ble |l] These parameters have been calculated [5] based 
on the 6-oscillator model fitted to the tabulated optical 
data given in Ref. [23, . 



j 


a;, [eV] 


9j [eV] 




1 


3.05 


0.75 


7.091 


2 


4.15 


1.85 


41.46 


3 


5.4 


1.0 


2.700 


4 


8.5 


7.0 


154.7 


5 


13.5 


6.0 


44.55 


6 


21.5 


9.0 


309.6 



sen the values Xt/R = 1.52, 5.94 since they correspond, 
e.g., to a sphere of radius R — 5/im, which is small but 
still relevant to experiments. The temperature is cho- 
sen as r = 300° i^: and T = 77°K, yielding At = 7.6^m 
and At = 29.7//m, respectively. These two temperatures, 
corresponding to room temperature and the boiling point 
of molecular nitrogen N2 respectively, can readily be ac- 
cessed in experiments. 

Below, we employ more detailed models for the mate- 
rial response to calculate the Casimir energies at higher 
temperatures. More specifically, we consider generalized 
plasma and Drude models, which take into account the 
interband transitions of core electrons that are described 
by a set of oscillators with nonzero resonant frequencies. 
The generalized plasma model has the dielectric permit- 
tivity 



ep{icK) = 1 + ( ) + ec{icK) 



and the generalized Drude is described by 



ep{icK) = 1 + 



(Apfc)^ -|- TTCK/ap 



with 



K 



(32) 



(33) 



(34) 



Here K is the number of oscillators, jj are the oscillator 
strengths, gj are the relaxation frequency and lOj ^ 



TABLE I: Oscillator parameters for gold. Calculated in Ref. 
t5j by fitting 6 oscillators to tabulated optical data [23] . 

We first calculated the n = term of the sum of 
Eq. ([1]) analytically, based on the expressions given in 
Section IIVI and then calculate the other terms numeri- 
cally as explained previously. This allows us to calculate 
the Casimir free energy for very short separations. As 
explained above, large values of imax should be used at 
short separations but in this limit the numerical evalua- 
tion of the determinant in Eq. ([I]) is cumbersome. This 
problem disappears as the temperature in increased be- 
cause then the second Matsubara wave vector ki becomes 
larger and thus no calculation for too small values of 
K > will be necessary. For the purpose of this paper, 
we calculate the relevant energies for distances as short as 
R/d — 0.95 at the two different temperatures using the 
generalized form of the Drude and the plasma model, see 
Eqs. ([32|, ([33f. 

Figr]4 shows the ratio of the Casimir free energy to 
the energy at T = for the generalized plasma model 
and for the generalized Drude model at T = 300° iiT and 
T = 77° K. While the Casimir free energy at T = 300° 
is always larger than at T = 77° K for the plasma model, 
we find that for the Drude model the Casimir free en- 
ergies at T = 77° K and T = 300° if cross each other 
around R/d — 0.7. For R/d> 0.7, the Casimir free en- 
ergy corresponding to T = 77° K becomes larger than 
the one for T = 300° if. For the sphere plate geome- 
try, indications of negative entropy have recently been 
reported [T7]. The ratio shown in Fig. |4] can be ex- 
pressed as E{T)/E{0) = 1 - TS/E{0) where S is the 
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FIG. 4; Ratio of the Casimir free energy to the energy at 
T = against R/d for the generalized Drude model (open 
symbols) and the generalized plasma model (filled symbols) 
at T = 300° if (squares) and at T = 77° K (circles) for a 
sphere of radius R = 5/^m. 
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FIG. 5: The ratio of Casimir free energies for the generalized 
plasma and generalized Drude model at T = 300° if (circles), 
T = 77°K (triangles) and T = (squares). 



VI. CONCLUSION 



entropy associated with the Casimir free energy E{T). 
Hence, a ratio E{T) / E{0) < 1 implies a negative en- 
tropy since E{0) < 0. Our results clearly show that for 
the Drude model the entropy indeed becomes negative for 
sufficiently small separations. However, for the plasma 
model our results of Fig. [4] indicate a positive entropy. 

Above we showed that at large separations the ratio 
of the Casimir energy for the plasma and the Drude 
model varies between 3/2 (for small Xp/R) and 1 (for 
large Xp/R) for zero and finite temperatures. At shorter 
separations the ratio is expected to tend to one since 
at high frequencies the plasma and Drude model be- 
come identical. It is interesting to observe how the ra- 
tio tends to one with decreasing separation as a func- 
tion of temperature. Fig. [5] shows the Casimir free en- 
ergy for the plasma model divided by that for the Drude 
model at T = 300° T ^ 77° AT and T = 0°K. Since 
Xp/R = 0.0274 is small compared to one, the ratio tends 
to almost 3/2 at large separations. As shown in the fig- 
ure, with decreasing separation the ratio drops towards 
one very fast for T = 0. However, for T = 300° AT the 
ratio is larger than 3/2 for R/d < 0.7, goes through a 
maximum around R/d = 0.6 and finally starts dropping 
to one. The curve for T — 77° K also displays a slight 
maximum close to R/d — 0.15. A similar behavior with 
an extrema has been observed also in Ref. \17\ for a suffi- 
ciently large sphere at T = 300° A". The maxima occur at 
a distance that approximately corresponds to the ther- 
mal wavelength At with R/Xt = 0.66 and R/Xt = 0.17 
for T = 300° a: and T = 77° K, respectively. Since ther- 
mal photons of wavelength At mostly contribute to the 
energy at a separation d « At, the position of the max- 
ima suggests that thermal effects less strongly enhance 
the Drude energy than the plasma energy, presumably 
due to dissipation. 



We have shown in detail how the scattering approach 
for Casimir interactions can be applied to study corre- 
lations between effects due to geometry, material prop- 
erties and finite temperature. The experimentally most 
relevant geometry of a sphere and plate reveals interest- 
ing properties of the Casimir interaction that are absent 
for parallel plates and hence in the proximity force ap- 
proximation. These findings demonstrate an interplay 
between material properties and the finite size of the 
sphere. Our main results are as follows. At large sep- 
arations we observe both at zero and finite temperatures 
for the amplitude of the leading term of the energy dif- 
ferent results for perfect reflectors, Drude and plasma 
metals. The plasma model yields a non-universal ampli- 
tude that depends on the ratio of plasma wavelength to 
sphere radius. For the perfect reflector and Drude model 
the amplitudes are universal but for the latter it is re- 
duced by a factor of 2/3. This result is distinct from the 
interaction of two parallel plates, which at zero tempera- 
ture is asymptotically identical for the three material de- 
scriptions. The identification of the plasma wavelength 
with the penetration depth of a London superconduc- 
tor explains why the plasma model yields the asymptotic 
interaction for perfect reflectors and a Drude metals as 
limiting cases. 

Our numerical computations of the energy at smaller 
separations demonstrate further important differences 
between the plasma and Drude model and generaliza- 
tions thereof. We observe full agreement of the numeri- 
cal results with the asymptotic expansion at large sepa- 
rations that, however, is limited for the Drude model to 
extremely large distances. Hence, we conclude that for 
the Drude model higher order multipoles are more impor- 
tant than for the plasma model. At small separations the 
observed dependence of the difference between the exact 
and the PFA energies on the plasma wavelength demon- 
strates that geometry and material effects are correlated. 
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Our results at finite temperatures show that the Casimir 
energy for Drude metals changes non-monotonically with 
temperature, leading to a larger energy at T = 77°K than 
at T = 300°K at sufficiently small separations. We ob- 
serve a negative entropy associated with the Casimir free 
energy for the Drude model over a range of distances. 
This range increases when the temperature is decreased. 
Both non-monotonic temperature dependence and neg- 
ative entropy are not observed for the plasma model in 
the range of studied parameters. At finite temperatures, 
we find that the Casimir free energy for plasma metals 
is approximately 3/2 times the energy for Drude metals 
for separations d ^ Xt- 
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